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The X-ray edge problem of graphene with the Dirac fermion spectrum is studied. At half-filling 
the linear density of states suppresses the singular response of the Fermi liquid, while away from 
half-filling the singular features of the Fermi liquid reappear. The crossover behavior as a function 
of the Fermi energy is examined in detail. The exponent of the power-law absorption rate depends 
both on the intra- and inter-valley scattering, and it changes as a function of the Fermi energy, 
which may be tested experimentally. 
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I. INTRODUCTION 

Recently the graphenes are being studied intensively 
mainly due to its novel transport properties such as 
quantum Hall effect Many of such novel proper- 

ties are attributed to the Dirac fermion band structure 
of a graphene^ The density of states (Dos) is linear in 
energy, and vanishes at the Fermi energy at half-filling. 
Therefore, many properties of graphene at half-filling are 
expected to be markedly different from those of the or- 
dinary Fermi liquid (FL) which has a finite Dos at the 
Fermi energy. However, away from half-filling, the Dos of 
graphene does not vanish at the Fermi energy although it 
can be very small, and some features of FL are expected 
to emerge. 




FIG. 1: The dependence of power-law exponent a of graphene 
on the coupling constant g — pVo for N c — 8 (see text for 
definitions). From top to bottom, ^ = 0,0.4,0.6,1. Note 
that the dependence on g implies the dependence on the the 
Fermi energy. The dashed curve is for the Fermi liquid. 



The problem of the transient optical response of FL has 
a long and interesting history£&£ ' 10 i 11 The high-energy 
X-ray creates a deep core hole, and as a result the Fermi 
sea of conduction electrons reacts to a time- dependent 



localized potential. In terms of the optical response the 
conduction electrons react in two distinctive ways: the 
excitonic processes^ and the orthogonality catastrophe 
(OC)2. For an attractive potential stemming from the 
created hole, the excitonic processes tend to enhance the 
optical response, while OC effects suppress the optical 
response for either repulsive or attractive potential. It 
is crucial to include these two effects on equal footing. 
For Fermi liquids with a short-range core hole potential 
V(r) = —V5(r), the absorption intensity which is asymp- 
totically exact near threshold takes the following form. 10 



I(oj) ~ 6(u> — Wy) p 



- UJn 



E r . 



(1) 



where u)^ is the threshold energy which is the energy 
gap between the core level and the Fermi energy of FL 
with corrections from particlc-holc interaction. 9{lo) is 
the step function. The power-law exponent is given by^ 



S 



+ N s (~y, 5 = tan- 1 (irpV). (2) 



The constants p, E c , N s , S are the Dos at the Fermi 
energy, a high-energy cutoff, the spin degeneracy 2, and 
the scattering phase shift, respectively. 

The X-ray edge problem and the OC have the univer- 
sal character of the dynamic quantum impurity in gapless 
host system, and have been studied in diverse physical 
systems. The tunneling junction with a localized impu- 
rity having multiple scattering channel has been studied 
by Matveev and Larkii>i2. using an analogy with the X- 
ray edge problem. The X-ray edge problem and the OC 
in one-dimensional Luttinger liquid have been studied 
extensively in Ref . p^ll^fl^lT7lfl8j . Khveshchenko and 
Andersor*i2 studied the X-ray edge problem of the to- 
mographic Luttinger liquid which has a vanishing Dos 
at Fermi energy. Patton and Gellcr 20 have obtained the 
anomalous tunneling Dos of strongly correlated electron 
systems by mapping the original Hamiltonian to the one 
to which the method of X-ray edge problem can be ap- 
plied. We also note that the Kondo physics belongs to 



2 



the same class as the X-ray edge problem. A subclass of 
the Kondo physics problems where the conduction elec- 
trons have power law Dos (soft gap) at Fermi energy has 
been studied very extensively mainly in the context of 
the anisotropic superconductors and the dynamical mean 
field theory^^ 

At a sufficiently large filling graphene is expected to ex- 
hibit FL characters. On the other hand, at and near half- 
filling the conventional approaches to the optical response 
of FL cannot be applied directly due to the qualita- 
tive differences in band structure between graphene and 
FL. For example, the bosonization metho d 11 ! 24 cannot 
be applied because the local Green function of graphene 
cannot be represented by that of one-dimensional chiral 
fermions.— Another new feature of graphene, not present 
in the ordinary Fermi liquid, is the inter-valley scatter- 
ing. Here we employ the method of functional integrals 
to overcome these problems, and integrate out all con- 
duction electron degrees of freedom except for the one 
at the position of the immobile hole. As far as the local 
response is concerned, the resulting local action contains 
all the necessary information. 25,26 The analysis of the lo- 
cal action based on path integral approach enables us to 
study the FL and the graphene in a parallel way, high- 
lighting the similarities and the differences at the same 
time. 

In this paper we report theoretical results on the X-ray 
edge problem of graphene based on the massless Dirac 
fermion model. We find that the X-ray absorption inten- 
sity has the following power-law dependence (u> > uj%,) 



/OJ — UJ^\ a E c /UJ ~ UJ^\ 

V E n J + ~\ E c J 



'uj - uj* \i+N c (s/^y 



with the exponent 



a = -2— + N c (S/w) 

7T 



(3) 



(4) 



where 6± = tan -1 [7175(1/0 ± Vj.)] with the short-range 
intra- and inter- valley scattering strengths Vq and V\ 1 
respectively. In the basis of Dirac fermions (see below) , 
the scattering potential takes the form of the following 
4x4 matrix which is consistent with the symmetry of 
graphene i2£ 



V scat (r) = 6(f) (v h + Via x ® 7- B ) = 8{r)V. 



(5) 



The combination Vq ± V± in the phase shift 8± is nothing 
but the eigenvalue of the scattering matrix V, and the 
Green functions are also naturally decomposed with the 
corresponding eigenvectors, see Eq. (|22[) . The exponent a 
of the first term in Eq.([3]) can be negative, which implies 
that the absorption intensity is singular at the thresh- 
old energy uj, On the other hand, the second term is 
non-singular since the corresponding exponent is always 
positive. However, near half filling the second term of 
Eq.Q dominates since p oc fx, while the first term dom- 
inates away from half filling (the Fermi energy p is zero 



at half-filling) . The exponent a contains N c — 4 (8 if 
spin included) in contrast with the Fermi liquid case of 
N s = 2. Since the phase shifts S± depends on p, the 
exponent a also depends on the Fermi energy p, see 
FigfT] In graphene the electron density can be varied 
up to 10 13 cm~ 2 by controlling the gate potentials, 1,2,3 
and this opens up the possibility of the systematic mea- 
surement of the dependence of the exponent a on the 
electron density. As we discuss in Sec. VII our result for 
the X-ray absorption intensity has several other interest- 
ing experimental consequences. 

This paper is organized as follows: In Sec. II we de- 
scribe our model Hamiltonian. In Sec. Ill its local ac- 
tion is formulated within path integral formalism. The 
Green functions of Dirac fermion and the core hole are 
derived in Sees. IV and V. Using these Green functions 
we derive the X-ray absorption intensity in Sec. VI. The 
discussions are given in Sec. VII. 



II. MODEL HAMILTONIAN FOR GRAPHENE 

Graphene can be described by the following tight- 
binding Hamiltonian of 7r-bandSf22, 



Ho = -tJ2 [c* l (R)c B (R + r i )+h.c. 
R.,ri 



(6) 



where A, B denote sublattices, and R (with r^) spans 
a honeycomb lattice. The hopping parameter is t — 
— 3.033 eV, and the basis vectors are given by ai — 
£{^,-\), a 2 =£(0,l). The lattice spacing is 1 = 2.46A. 
Each site of A sublattice is connected to three B sublat- 
tice sites: = £(^,1), r 2 = £(^,-1), r 3 = 

t{— ttjjO). The energy eigenvalue is given by E±(k) = 

±*|A(k)|, where A(k) = e lkri + e ik ' r2 + e ik ' r3 . Two 
bands touch at 6 points located along the corners of 
Brillouin zone, and among them only two are distinct: 
K+ = ^ (i, j^=) and K_ = -K+. These two points 

will be denoted as K and K', respectively. From now on 
we will set the lattice constant £ = 1 and h = 1. 

At low energy the above Hamiltonian reduces to the ef- 
fective Hamiltonian of the massless Dirac fermion, which 
is valid within a few eV of the K and K' points 5 -* 6 . In a 
4-component notation involving two from sublattice A, B 
and two from valley indices K±, the Hamiltonian has the 
fori 



,29 



Hn = d 



■.^{x)(v{-id x )a x +v(-id y )a y ^(x), (7) 



where v oc t is the velocity of Dirac fermion and ^ = 
[ca+, cb+, ca-, cb~] T is the 4-component spinor (spin 
indices suppressed). The gamma matrices are a x — 



■h 



v" = 



o~y (g> (— t 3 ), j3 — a z <S> I2, and they sat- 
isfy the standard gamma matrix algebra {a 1 , a^} = 25 ZJ . 
Tj, <7j are two independent copies of Pauli matrices acting 
on the sublattice and the valley space, respectively. 
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The Matsubara Green function of Dirac fermion with 
Fermi energy fi is given by 

G,6(ie,k) = -(¥„(»£, k)*J(»e,k)) 

{ 1 \ 



ie + fi — V J2j =x ,y kjai J ab 
ie + v k J a J \ 



(ie + fi) 2 — v 2 \i 2 



(8) 



ab 



where a,b are the collective sublattice/valley indices, ie 
is the imaginary frequency, and k J is the j-th component 
of wavevector k. In the last line of the above equation 
the gamma matrix algebra has been used. 

In our model the core potential is taken to be isotropic 
since the core hole belongs to the s-level. Also, it is as- 
sumed that the core potential is screened instantly, and 
this assumption should be valid^ near the threshold en- 
ergy ui T . Thus the deep core hole can be described by 



H; 



dU^(x = 0)V^(x = 0), 



(9) 



with the scattering potential V given by Eq.©. d is the 
annihilation operator of the deep core electron and its 
Hamiltonian is given by Ed d^d, where Ed < is the core 
level energy. The (bare) Matsubara Green function of 
the core hole is given by 



D (ie) 



1 



ie + ii- Ed 
The interaction with X-ray is described by 
H x = J2 M Q>k *t (kjde"*" 4 + H.C. 



(10) 



(11) 



M a ,]i is the X-ray matrix element. Usually the matrix el- 
ement is very weakly dependent on k and its dependence 
will be ignored. The dependence on the index a will be 
also ignored. In this paper the interaction between Dirac 
fcrmions is not taken into account. 



III. FORMULATION IN TERMS OF LOCAL 
ACTION 

The local Matsubara Green function of Dirac fermions 
of graphene with the Fermi energy fj, is given by 

G ab (r - t') = -(*„(£ = 0, r) ^l(x = 0, r')) 
-A 1 M 2k e- i <—'\ie + l i) 

~ ° ab 2^ J ( 27r )2 ( ie + M )2 - V 2 k 2 ■ ^ 

Note that in the last line of the above equation, the term 
depending on momentum linearly [ a J k J , see Eq.® ] 



vanished upon wavevector integration. If |r — t'\ ^> r c , 
where r c is a short time cutoff, Eq. lfT!?]) becomes at zero 
temperature 



G a b(l~ - t') = S ab G = 5 a b(GFL + Gt>), 

1 1 



G d (t-t') 
Gfl(t-t') = — 



t' — t) er^'^^ 



4irv 2 (r - r') 2 Attv 2 (t - r') 2 ' 
^ (13) 



where 



Attv 2 ' 



(14) 



which plays the role of Dos at the Fermi energy. Note 
that it is proportional to fi. Thus the local Green func- 
tion consists of the Fermi liquid component Gfl and the 
pure Dirac component Gd- In the long time/low energy 
limit the Dirac component is much less singular than the 
FL component. In fact, the direct Feynman diagram cal- 
culations at half-filling show the logarithmic singularities 
originating from Gfl do not exist, so that the perturba- 
tion theory converges well at low energy. For details, see 
Sec. VII. 

To obtain the action of the local degrees of freedom 
we have to integrate out the bulk degrees of freedom ex- 
cept for the one at x = 0. For this step, we use the 
method of Dirac delta function in Grassman variable in 
the following form 

r 4 

1 = / D[ Va ,f] a } Y[ 6(n a - ® a ,x=a)S(fja - K, x =o) 

J a=l 

= / A a ,A a ] exp [J^ \ a (T)(7] a (T) -* o (0,t)) 

+ (Va(r)-¥ a (0,T))X a (r)], 

(15) 

where A a is a Grassman Lagrangian multiplier. The local 
action of the local degrees of freedom r] a (T), fj a (T) can be 
obtained by integrating out \1/ Q and A a successively. 

The resulting local action of r\ a and d in imaginary 
time is 

S[r),d] =-E/ dTdT'f) a {T)G- b HT-T>) m {T') 

a,b 

+ J dr d(T)d(T)fj(T)Vr)(T) + J drd{d T -uj T )d, (16) 

where lot = M — Ed is the (bare) threshold energy and 
G~ b is the inverse of the local Green function Eq. (fT!2|) of 
the Dirac fermion. The local action Eq. ([To]) is of a very 
general form, so that it can be applied to other cases with 
different band structures, including FL. The specific band 
structure is reflected only in the local Green function 

G ab (T-T>). 
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IV. THE GREEN FUNCTION OF DIRAC 
FERMION 



are proportional to Gfl, while g Q J 1 do not. The explicit 
forms of these new Green functions are 



In the presence of a following time- dependent potential 

V(t) = -V9(n - t)9(t - t 2 ), (17) 

the local Matsubara Green function of Dirac fermions is 
given by (£,£' are imaginary times) 



ff(c-0 = <ci 



i 



G~ 



V 



10, 



(18) 



which is a 4x4 matrix Green function, and this Green 
function encapsulates the excitonic contribution to the 
X-ray response function. From the above definition of 
g one can show that it satisfies the following (matrix) 
integral equation. 

9(t,0 = G(t;-e)+ P dTG(S,-T)(-V)g(T,a, (19) 



FL 1 (-p) cos 2 5+ (t x -t 2 
5o,i( r i ~ T 2) ~ o — Z Z — 

2 71 - T 2 



1 {-p) cos 2 5- (t\- t 2 
2 



9o(n r 2 ) 

ff?(n - t- 2 ) 



n - r 2 

1 1 

4ttw 2 (n — r 2 ) 



47TW 2 



, Vi 1 

47ri> 2 t c 



1 



(ri - r 2 ) 2 ' 



(25) 



where p is to be understood as the renormalized Dos at 
the Fermi energy. From these it follows that 



5o,i « 5o,i + 5o,i- 



(26) 



V. 



CORE-HOLE GREEN FUNCTIONS 



where G is the bare local Green function, Ea.(|13p. which 
is diagonal. To solve Eq. (fl"9|) we choose an ansatz 



(20) 



g = gah + gio- x ® t x , 



and take the combinations of go ± 31 ■ Then it can be 
shown that each of go ± g\ satisfies the following equation 

(g ±g 1 )=G+(-V TV 1 ) P drG((-r)(g ± gi ). (21) 

J T 2 

When G is of the FL form Gfl(t) = —p/r, the solution 
of Eq. (|21[) which is asymptotically exact in the long time 
limit was obtained by Nozieres and De Dominicis (ND) 10 : 



( 50 ±3i) ND (£,O = [cos 2 5 ± ]G FL (£-£' 

-r 2 ) 

with the scattering phase shift given by 
6± =tan" 1 [7r(V r ±V r i)p]. 



(22) 



(23) 



This result cannot be obtained perturbatively. In view 
of the fact G = Gfl + Gd, the Green functions go.i can 
be also written as a sum of 



5o,i 



ND 

5o,i 



Sgo,i- 



(24) 



where <5go,i can be computed perturbatively since Gd is 
much less singular than Gfl in the long time limit. Us- 
ing the above form of go.i in Eq. (|21|) and carrying out the 



integrals carefully we find that Sgo,i ~ ^o.iSoT + A),iGd 
where ^4o,i are small multiplicative renormalization con- 
stants of the order of Vo,i/t and Bq i are also multiplica- 
tive renormalization constants. We define, for conve- 
nience, new Green functions g^\ — g^Y + ^a,i9oY^ ana ^ 
g J 1 = £?o,iGd- Note g$\ vanish as p — > since they 



The photoemission spectra is related to the core hole 
Green function which captures the physics of OC. 



D(t u t 2 ) = (T T d(n)d(T 2 )). 



(27) 



At temperature much lower than lot this Green function 
is non- vanishing only for t± > t 2 . Carrying out the path 
integral over d in Eq. (|16p we obtain 



D{t u t 2 ) 



JD[v}< 



^ T2 



In) 



fDM> 



(28) 



\ti) is Green function (1st order in time) in 



where S v is the first term of Eq. (Tl7a|) . In fact, a term 
-j(5(t— t') should have been included in S^. It stems from 
the determinant of c?-path integral, but it turns out not 
to affect the long-time behavior, so has been dropped. 

the presence of external potential f)Vi], which is easily 
calculated to be exp[— dT r fj(— V)rf\. Now the path 
integral over 77 can be done. Recalling G = Gfl + Gd, 
the result of the path integral can be expressed as 

D(n,T 2 ) = e THn(l-G FL y) e Trln(l- I ^)_ 



The first factor of Eq. ([29)1 is the core hole Green function 
of the Fermi liquid component, and this can be calculated 
exactly by using the ND solution Eq . (|2"2"]) The sec- 
ond factor does not generate any infrared singularities 
and it can be evaluated by a perturbative expansion in 
V. Expanding the second factor up to the second order 
in V we obtain 



D{T U T 2 ) ~ e -"r(«-7* 



1 



x exp 



N c (vj+vr 

4 



(4ttw 2 ) 2 ) (n - T 2 ) 2 



,(30) 
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where is the renormalized threshold energy-Z. and <5 
is the exponent given in Eq.(|3]). N c — 4 is the num- 
ber of the species of Dirac fermions, and its presence 
originates from the trace over spinor space. The result 
Eq. (f30|) depends critically on band filling. At half-filling 
we have 6 = 0, and the correction from the scattering 
is negligible in the long time limit. But away from the 
half- filling the power law dependence dominates. Eq. (f3"0)) 
also implies that the overlap between ground state wave- 
functions with and without the core hole is given by 



1(01*01 



g+jJVcd-) 2 ln(w min T e )-const. (u, 



(31) 



where w m ; n is the nv/ L with the system size L. At half- 
filling the first term in the exponent, which is of the FL 
type, vanishes since 6 = 0. However, the second term 
is finite and the overlap becomes almost independent of 
the system size. This is is consistent with the result of 
Ref. [3l[ which was obtained numerically. 

The hole spectral function can be obtained from (£ is 
an infinitesimally small positive number) 



A h {u) 



-Im 



dre iuT D{T) 



In the long time limit (or in the vicinity of id = wj) 

\2i I T c c l 



(32) 



A h (uj) ~ 9(id - Ut) sm[nN c (S/TrY 



[(w-^jrj 1 -^*)" 



C 2 T, 



(33) 



where (r(x) is the Gamma function) 
Cl = T(l - N c (6/ir) 2 ), 



C2 = 



V 2 



^r(i-jy e (W) 

(4ttw 2 ) 2 J (1 + N c (8/%) 2 )N c (6/%) 



(34) 



The above expressions of are valid as long as 1 > 
N c (6/n) 2 . If N c {6/ir) 2 > i, the r-integral of Eq.® 
diverges in the short time limit r — > 0. However, we have 
to remember that the expression Eq. (|30"|) of D(t) is valid 
asymptotically only in the long time limit. Therefore, 
this divergence is not physically relevant, and it can be 
treated by introducing an explicit short time cutoff in 
the integral or by an analytic continuation of Gamma 
functions appearing in c\ and C2. 

At half-filling we have 6 = 0, and the core-hole Green 
function Eq. (|3"D|) becomes 



D(n,T 2 ) 



-(ti— t 2 ) 



exp 



Nc 
' 4 



V 2 



(Anv 2 ) 2 ) (n - r 2 ) 2 
(35) 

Expanding the exponent of Eq. (|35p in the long time limit 
we obtain 



Ah,s=o{Lo) - tt6{lo - 4) 

— 9 (id — Wj)(w — U>t) 



v? 



(4wv 2 ) 2 



(36) 



Comparing with the result of Eq. (|33|) (6 ^ 0), we find 
that the anomalous power-law dependence has disap- 
peared. 



VI. X-RAY RESPONSE FUNCTION 

This is given by 

F(n,T 2 ) = ^2(T T d(n)r] a (Ti)fj b (T2)d(T2)), (37) 

a . b 

which is also non- vanishing only for t\ > t 2 . To compute 
this response function it is convenient to introduce the 

external source J a for r\ a , Sj = - J dr(^J2 a JaVa+fjaJa^j , 
into the local action. Carrying out d-integration and J- 
differentiation we arrive at 



a, 6 J L J 



J D[rj\ e~ s W~ Sj ~f dr vVn 



j — o 



(-1) 



DetfG- 1 - V] 
Det[G" 



i ^9ab(Tl ~ T 2 ), 

' „ h 



(38) 



where g a b(T\ — t 2 ) is the (matrix) local Green function 
of the Dirac fcrmion of Eq. (fT8"]) . We also note that 



Det[G~-V] • 



Dct[G _1 ] 



is the core hole Green function of Eq. 



The path integral formulation automatically yields that 
the excitonic effect (represented by g a b) and the OC (rep- 
resented by the core hole Green function) contribute in 
multiplicative way. Therefore, we obtain 

F(n-T 2 ) = -6(n-T 2 )D(n-r 2 ) Y.9ab(Ti-T 2 ). (39) 

a . b 

From Eq.® it follows that 

^2gab(n - t 2 ) = 4(g + gi) 

= 4kg FL +.9f L ) + (.9o D +5?)l-(40) 



a , b 



Now the absorption intensity at T = can be computed 
from 



I(w) 



1. 



-Im 



dTe iuJT F(r) 



iu—^LoA-i^ 



(41) 



Using Eqs. (|25l3QI39[) we find for the absorption rate 

/ = /fl+/d, (42) 

where 

- 2 ^±+N c (S/tt) 2 



Jpl(w) ~ C?FL 9 (id - id?) P 



E, 



Iu(id) ~ C?D 9 (id — id?) 



pE c I id — id^ 



(1 



l + 7V c (<5/7r) 2 



, (43) 
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with the coefficients 

c?fl = cos 2 6+ T(2^± - N c (6/ir) 2 ) sin [tt(2^ - ^(^/tt) 2 )] , 



function of graphene (at half-filling and at zero temper- 
ature) is given by 



d D = sin [N c {6/n) 2 ] 



r(l - N c (S/tt) 2 ) 
[l+N c {8/Tif]N c {5/nf 



(44) 



The above expressions of c?fl and do are valid as long as 
1 > N c (S/tt) 2 . In the case of 1 < N c (S/tt) 2 , d F L,D are 
given by an analytic continuation of Gamma functions as 
in the case of the core-hole spectral function. In particu- 
lar, c?fl is positive even if 2^- — N c {^) 2 becomes negative, 
owing to the identity 1 = T(l — x)T(x) sin(7rx)/7r. Note 
that only the exponent S + appears in Eq. (|42]) because 
the second term of 5o^( T i ~ T 2) which has a 6- exponent 
[see Eq. ([25|[40|) ] is canceled in the summation over chan- 
nels a, b of Ea. ([551) . In the above we have assumed that 
the X-ray matrix element is very weakly dependent on 
the sublattice and valley structure. Let us compare our 
result of the exponent of Eq.|g2J) with that of Ref.[H[. 
Eq.(9) of Ref. [13| is the leading term of the exact result 
Eq.(72) of Ref. [lOj. In our case only s-wave orbital scat- 
tering is relevant, and the effect of channels other than 
the orbital such as valley and sublattice is reflected in the 
factor N c . Therefore, our result is fully consistent with 
that of Ref. [3. 

The /fl(w) vanishes as half- filling is approached since 
p and 8± — > at the same time, while the 7d part re- 
mains finite. The exponent —2^- + N c (^) 2 is a sum of 
the exponents for the excitonic and OC processes, —2^- 
and N c (^) 2 , respectively Away from half-filling, if this 
exponent is negative, /fl is the dominant term since it 
diverges as u> — > lot, manifesting the X-ray singularity. 
The Id is finite and does not exhibit a singular behavior. 

At half-filling, even in the absence of the local poten- 
tial, the intensity is suppressed by a factor of lo — lo t 
near threshold, and this is a consequence of the effect of 
the linear Dos. The linear Dos completely suppresses the 
excitonic processes near the threshold, and the perturba- 
tive treatment is sufficient, as can be checked explicitly 
by a direct calculation of Feynman diagrams (see below) . 



VII. DISCUSSIONS 

The singular nature of X-ray edge problem of FL can 
be traced back to the logarithmic divergence of 1-loop 
particle-hole polarization function.- Due to these log- 
arithmic divergences the infinitely many Feynman dia- 
grams should be summed even in the weak coupling re- 
gion, and the summation leads to the anomalous power 
law dependence of response functions. Thus to under- 
stand the non-singular nature of response functions of 
graphene at half filling (/i = 0, p = 0), it suffices to cal- 
culate the 1-loop particle-hole polarization function and 
to show that it is non-singular. The 1-loop polarization 



-N r 



de f d 2 k 



2n J (27r) 2 (ie) 2 — t> 2 k 2 ie — iui + lot 



(45) 



where G(ie,k) is the matrix Green function of Eq.([8]) 
at fi = 0, and Do(ie — iui) is the bare core hole Green 
function, given in Eq. (fT0|) . In the last line of the above 
equation the term linear in k vanishes upon k integration. 
The integral is straightforward, and we obtain 



n<°>M 



1 



2ttw 2 



(ioo — lot) In 



( Wt 



E r 



, (46) 



where E c is an energy cutoff of the order t. As demon- 
strated in the above (after analytic continuation iio — > 
lo + i£), the logarithmic singularity near the threshold 
\vl(u)t — is completely suppressed by the prefactor 
lj — lot which stems from the linearly vanishing density of 
states of graphene. In the case of FL, the prefactor lo—lot 
is replaced by the (finite) density of states p. Therefore, 
for graphene at half-filling, the higher order corrections 
are not necessary if we are only interested in the vicinity 
of threshold. 

It would be interesting to observe in X-ray experiments 
of a single layer graphen o 32 ' 33 how the exponent a of 
Eq.Q depends on the electron density. In graphene the 
carrier density can be varied up to 10 13 cm -2 by changing 
external gates^ 2 * 3 -. This corresponds to a change in the 
Fermi energy up to a few eV. The typical X-ray energy 
range is given by the energy difference between the core 
electron in the Is level and Dirac fermions residing in the 
7T band, which is in the region of 280 ~ 290 eV2&2£ One 
of the key experimental features would be that the ab- 
sorption spectrum diverges at the threshold energy lo t , 
which would correspond to a negative exponent a in the 
first term of Eq.Q (the second term of Eq.([3]) is zero at 
the at the threshold energy lo = lo t ). In the Fermi liquid 
the exponent decreases monotonically with the electron 
density, while in graphene it decreases initially and then 
increases and becomes positive, see Fig.l. When the ex- 
ponent becomes positive the absorption spectrum does 
not diverge at the threshold. Experimental verification 
of this change of sign of the exponent in graphene would 
be intersting. The experimental measurement of a would 
also provide useful information about the strength of in- 
tra and inter valley scattering strengths Vq and V\. 

In this paper, the hole was assumed to be immobile 
(or of infinite mass). 20 ' 34 ' 35 ' 36 The finite mass effect and 
the finite temperature 3 ^ are expected to smear out the 
singular behavior to some extent, and these factors are to 
be taken into account for the detailed comparison with 
experimental data. 
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